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Abstract 

p I | Solutions to partial differential equations with highly oscillatory, large random potential 

have been shown to converge either to homogenized, deterministic limits or to stochastic limits 
depending on the statistical properties of the potential. In this paper we consider a large class of 
piece- wise constant potentials and precisely describe how the limit depends on the the correlation 
properties of the potential and on spatial dimension d > 3. The derivations are based on a 
Feynman-Kac probabilistic representation and on an invariance principle for Brownian motion 

£SJ ' in a random scenery. 

>' 

in 

§ ■ 1 Introduction 

in 

Partial differential equations with small scale structures arise in many aspects of physics and applied 
science, and homogenization theory has proved to be useful both from the theoretical and numerical 
point of view. Depending on the statistical property of the small scale coefficients, the limit, as 
the size of the fluctuations tends to zero, could be either deterministic or random. Understanding 
the dependence of the limiting equation on the coefficients is therefore important in many practical 
settings. 

Parabolic equations (including heat Anderson or Schrodinger models) with large, random Gaus- 
sian potential are analyzed in [1, 2, 3, 19, 20] using Duhamel infinite series expansions and combina- 
torial techniques. The limiting equation is either a deterministic, homogenized, parabolic equation 
or a stochastic partial differential equation depending on the correlation properties of the (possi- 
bly time-dependent) potential; see [3]. In [16, 9], a one dimensional heat equation with singular 
time-dependent coefficients is studied by probabilistic (Feynman-Kac formula) and analytical tech- 
niques, respectively. The limiting equation is also deterministic or random depending on whether 
the coefficients have microstructure in both temporal and spatial variables. In this paper, we fol- 
low the Feynman-Kac approach to analyze a parabolic equation in dimension d > 3 of the form 
dtu e = Au e + iV e u e , where V e is a large, time- independent, random potential, which is not neces- 
sarily Gaussian. An imaginary potential is introduced to obtain a uniform bound on the energy 
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of the solution, which considerably helps in the analysis of exponential functionals of Brownian 
motion and of the limit as e — > 0. The corresponding heat equation (with iV £ replaced by V £ ) might 
be analyzed using techniques in [9] but is not considered here. The above scalar equation may be 
recast as the system 

in which V £ might model a conservative process of interaction between two populations and find 
applications in biological or social sciences. 

The Feynman-Kac representation is applied in [15] to a one-dimensional heat equation with time- 
independent coefficients. Using the same probabilistic representation, as well as auxiliary problems, 
[14] addresses homogenization of divergence-form operators with lower-order random potentials 
written in terms of gradients of stationary processes. In [10], the authors study the Feynman- 
Kac formula for a heat equation driven by a fractional white noise in all spatial dimensions. This 
corresponds to V £ with a sufficiently long-range correlation function (i.e., sufficiently slow decay). 
In [13], the heat equation with long-range correlated Gaussian potential is studied with a similar 
type of limiting equation as in [10]. 

If we consider the random integral appearing in the Feynman-Kac formula, i.e., L V £ (B s )ds, 
where V £ is the random potential, and B s is some independent Brownian motion, we obtain the 
so-called problem of Brownian motion in a random scenery. In the discrete case, this corresponds 
to the Kesten-Spitzer model of random walk in a random scenery. An invariance principle has been 
proved in [11, 6] in all spatial dimensions. A continuous version for a specific type of potential 
has been analyzed in [17], primarily for d = 1,2. In [12], the authors analyze the case where the 
Brownian motion is replaced by some reversible stationary Markov process and prove a central limit 
result. In the aforementioned examples, the random scenery is independent of the random walk. 
This is different from the so-called random walk in random environment, where the underlying 
environment indeed affects the random walk, e.g., in the random conductance and percolation 
clusters model. For Brownian motion in random scenery, in addition to central limit theorem, most 
of the existing results are of large deviation type, and we refer to [18] and references therein. 

In [7], we prove an invariance principle for L V £ (B s )ds in all dimensions when the potential is 
stationary Gaussian or Poissonian. Although we do not do it here, the limiting results obtained in 
this paper can be extended to those case without much effort. Instead, we adopt the specific class 
of potentials introduced in [17], which allows us to treat the short-range and long-range cases in a 
uniform way, and prove a weak convergence result when d > 3 and obtain the asymptotic limit of 
our PDE solution. As mentioned above, the limiting equation turns to be either deterministic or 
random depending on the correlation properties of the potential V £ . 

2 Problem setup 

We are interested in the parabolic equation 



d t u £ = ^Au £ + i^V(^)u £ (2.1) 



with initial condition u £ (0,x) = f(x), x £ R , d > 3. Here V is a stationary random field with 
either short- or long-range correlation, which will be defined later. We wish to analyze the limiting 
behavior of u £ as e — > and its dependence on the random field. It turns out that as long as 
the covariance function of the random field does not decay too slowly, then homogenization and 
convergence to a deterministic solution occurs, while when the random fields are sufficiently long- 
range, the limiting equation is a stochastic PDE. Correspondingly, the large, highly oscillatory 
random field -£-V(§) is either replaced by a constant or by a Gaussian noise, respectively, in the 
limiting equation. 

The random field we consider is of the form V(x) = & x +in, where {£&, k £ Z, d } is a sequence 
of stationary random variables, U is uniformly distributed in [0, l) d and independent of {£&}, and 
[x] := ([xi], ..., [xd]) with [xi\ denoting the largest integer smaller than or equal to X{. Defining 
r{k) = EtfoCfc}, we verify that R{x) := E{V(y)V(y + x)} = E fe K*0Ilii(l " l*i " h\)l\ Xl - ktl< i, 
and thus, the decay of r(k) as \k\ — > oo determines the decay of R(x) as \x\ — > oo. 

The three cases of £& we consider are the following: 

1. Short-range correlation: {£&} are i-i-d random variables with zero mean and unit variance, 
and for technical reasons we assume E{£(j} < oo. In this case, 

d 

R(x) = l[(l-\x l \)l lXil<1 . (2.2) 

i=l 

2. Long-range correlation: 

• Case 1: {£&} is a stationary Gaussian sequence with zero mean and unit variance. We 
assume its covariance function 

r(k) = E{Mfc} ~ ^ (2-3) 

for some a € (2, d) as \k\ — > oo. In this case, 

d 



R(x) =Y,r{k) J](l " fa ~ ki\)l\ Xi - ki \<i ~ r^ (2-4) 



as \x\ —7- oo. 



Case 2: {£& = $(%)}, where {r/^} is a stationary Gaussian sequence with zero mean 
and unit variance. We assume its covariance function p(k) = E{ryo??fc} satisfies \p(k)\ < 
CJ\i—i min(l, TjAa-) for some constant C and 

p{k) ~ n^ \ n (2 ' 5) 

lli=i l^r 

as minj = i d \ki\ —> oo for some on G (0, 1) and constant c^ > 0. The function $ satisfies 
E{$(r/ fc ) 2 } < oo, V = E{$(r/ fe )} = and V\ = E{<S>(r] k )r] k } / 0. We further assume that 
«:=£ii«*e(0,2). 



Since r(k) = E{£o£fc}> we will see that r{k) ~ V±p(k) as |fe| —> oo, so 



i?(x) = J>(fc) J](l " \xi ~ h\)l\ x ^ kl \<i ~ -g 1 , Cd | (2.6) 



fc i=l llj=ll x i| 

as mirij =1 ^ |xj| — > oo. 

Remark 2.1. In the short-range case, the assumption E{£q} < oo is used in proving the tightness 
results of a family of processes in C([0, T]), which we do not necessarily need in the proof of Theorem 
3.11. 

Remark 2.2. In both long-range correlation cases, we assume a < d, which implies that R(x) is not 
integrable. 

In the first two cases, we define the positive constant that will appear in the homogenized 
equation 

-(wjW'-(;^>jLhB*) J >« ™ 

with R(£) = f Rd R{x)e~^ x dx. 

The scaling constant (3 is defined as 

J 1, short-range correlation, long-range correlation case 1, 
(f, long-range correlation case 2. 

In the short-range correlation case, the random field V{x) is mixing in the following sense. For 
two Borel sets A,Bc Mr, we denote by J- a and J~b the sub-cr algebras generated by the field V(x) 
for x £ A and x € B, respectively. There exists a mixing coefficient cp(r) such that 

\Cor( V ,()\<<p(2d(A,B)) (2.9) 

for all square integrable random variables rj and £ that are J- a an d J~b measurable respectively. 
The multiplicative factor 2 is only here for convenience. In the short-range correlation case, when 
|x| is large enough, V{x + y) is independent from V(y), so the mixing coefficient <j>(r) can be chosen 
as a positive, decaying function with compact support in [0, oo). We will use this in the calculation 
of the fourth moment of V(x). 

We also note that since the potential V(x) is discontinuous and possibly unbounded, (2.1) is 
not always solvable in a classical sense. We show that the solution defined by the Feynman-Kac 
formula is a weak solution in the following proposition. 



Proposition 2.3. Consider the equation 



u t = -Au + iV(x)u (2.10) 



with initial condition u(0,x) = f(x) G C&(R ). Let us define 

u(t, x) = E B {f(x + B t ) exp(i I V(x + B s )ds)}, (2.11) 

Jo 
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where E# denotes the expectation with respect to the Brownian motion starting from the origin. 
Then for any <p G C£°(R ), we have 



t r -i rt 



u(t,x)ip(x)dx = I f(x)(p(x)dx+ I I u(s,x)-Aip(x)dxds+ I I u(s,x)V(x)ip(x)dxds (2.12) 
d JR d Jo JR d 2 Jo JR d 

almost surely, i.e., the Feynman-Kac solution u(t,x) is a weak solution to (2.10) almost surely. 
Proof. Recall V(x) = £,[ x +u]> fi x an Y S,M > 0, define 

Vs.m(x) = / 4>&{x-y)V(y)l\ y \ <M dy, (2.13) 

JR d 

where <j>$ is a family of compactly supported mollifier. For each realization, since V{y)l\ y \ < j^ is 
bounded, V$ : m is bounded and as long as V is continuous at x, we have Vs t M(x) —> V(x) as 
S — > 0, M — > oo. In addition, V$ y M is smooth, so for the equation dtUs^M = \^ U 5,M + K5,Af^<5,A/ 
with initial condition u$ t ]\i(0,x) = f(x), we have its classical solution given by the Feynman-Kac 
formula 

u s , M (t, x) = E B {f(x + B t ) exp(i I V s ,m{x + B s )ds)}, (2.14) 

Jo 

and it is straightforward to check that u$ t M(t,x) — > u(t,x) as 5 — > 0,M — > oo by Lebesgue 
dominated convergence theorem. Since u^m is also the weak solution, we have 

us M(t,x)ip(x)dx = / f(x)(p(x)dx+ / / u$^m (s,x)-Aip(x)dxds 



t Jo jRd 2 (2.15) 

/ us,M(s,x)Vs : M(x)ip(x)dxds. 

JR d 
Let 5 — > 0, M — > oo, we finish the proof. □ 

Since V(x) may be unbounded, proving uniqueness of the solution to (2.10) is a difficult task. 
Such a task becomes easy when the equation is posed on a bounded domain since V is then 
bounded almost surely. But calculations with the corresponding Brownian motion on bounded 
domains involve standard complications which we wish to avoid here. When we refer to "the" 
solution to (2.1), we therefore mean the weak solution given by the Feynman-Kac probabilistic 
representation in the rest of the paper. 

In the following sections, we deal with the short-range correlation and long-range correlation 
cases in turn. Let us introduce some notation. We have randomness coming from two sources, the 
random potential V from the differential equation and the Brownian motion from the Feynman- 
Kac formula, and consider the usual product probability space. We use B = {Bt} to denote the 
Brownian motion starting from the origin and E# is used to denote the expectation only with 
respect to B. The expectation in the whole product probability space is denoted by E. iV(/x, a 2 ) is 
the Gaussian distribution with mean \i and variance a 2 , and qt{x) denotes the density function of 
N{0, t). We use a Ab = min(a, b) and a <b when there exists a constant C such that a < Cb. 



3 Short-range correlation 

Since we rely on the Feynman-Kac formula, the integral - L V(-f-)ds invites us to look at Brownian 
motion in a random scenary. Heuristically, V(x) is piecewise constant for each realization, so the 
integral is a sum of independent random variables, and we would expect a central limit theorem to 
hold. 

We first prove an invariance principle for Brownian motion in a random scenery. Using that 
weak convergence result together with the a priori bound for u e , we pass to the limit as e — > to 
derive the homogenized equation. 

3.1 Invariance principle for Brownian motion in a random scenery when d > 3 

Define X £ (t) = - J V(—)ds, the main result is the following theorem. 
Theorem 3.1. In the short-range correlation case, 

X £ (t) => a d Wt (3.1) 

weakly in C([0,T]), where Wt is a standard Brownian motion. 

To prove it, by scaling property of Brownian motion, we know X e (t) has the same distribution 
as e/n V{B s )ds, and defining 

p t {k) = \ l[B s +u]=kds, (3.2) 

Jo 

it remains to prove the weak convergence of X e (t) = e ^2 k £, k Pt/s 2 (&)■ This is divided into two steps. 

We first prove the convergence of finite dimensional distributions and then prove tightness. In the 

following, we do not distinguish between 1 f V(-f-)ds and £ J V(B s )ds. 

3.1.1 Finite dimensional distributions 

Let X(t) = f V(B s )ds = ^2k^kPt(k), we first prove two lemmas, which combined show that the 
variance of X e (t) conditioning on B, U converges in L 2 . 

Lemma 3.2. ^{Yl k Pt( k ) 2 } ~ cr d t - 
Lemma 3.3. E{(J2 k Pt(k) 2 ) 2 } ~ ajt 2 . 

Proof of Lemma 3.2. Since E{J2 k Pt(k) 2 } = K{X(t) 2 }, we have 

E{V Pl (if}=E{(/ V(B s )ds) 2 }= I I E{R(B S - B u )}dsdu 

u JO JO Jo 



1 /■*/■*. 

E{exp(i^(5 s - B u ))}R(£)d£dsdu 



(^) d Jo Jo 

: -i- J / / / exp(-ht\ 2 \s-u\)R(0dZdsdu 

(27r) c( Jo J J Rd 2 

1 L^wi'-w 11 --*^)* 



(3.3) 



(2vr) c 



Because d > 3, -R(£)/l£| 2 1S integrable, then by the Lebesgue dominated convergence theorem, we 
have 



^M^/.^/ t «* 



as t — > oo. □ 



Proof of Lemma 3.3. Since pt(k) is independent of the choice of £&, we can assume that £& 
have Gaussian distribution, in which case we have E{(^2 k p t (k) 2 ) 2 } = |E{X(£) 4 }, and the fourth 
moment can be calculated as 



f 4 
E{X(i) 4 } = 24 / E{[[V(B Si )}ds, (3.4) 

J St ._i 



' s * i=l 

where we denote 5f = {0 < si < ... < S4 < t}. It turns out that 

4 
E[]%) =P([si + ^] = [x 2 + U], [x 3 + 17] = [x 4 + t/]) 

+P([xi + 17] = [x 3 + U], [x 2 + U} = [x 4 + U]) 

+P([X1 + E/] = [X 4 + C/], [X 2 + C/] = [X 3 + U}), 



(3.5) 



so we write E{X(i) 4 } = 24((i) + (ii) + (w)), where 

(i) = / P([5 S1 + C/] = [J3 S2 + C7], [fi^ + U) = [B S4 + U})ds, (3.6) 

J St 

(ii) = / P([5 S1 + C/] = [5,3 + [/], [B S2 + 17] = [B S4 + U])ds, (3.7) 

(m) = f F([B S1 +U} = [B Si + U], [B S2 + U} = [B s , + U])ds. (3.8) 

J S t 

We first estimate (ii), (Hi). Note that P([x + U] = [y + U]) = R(x — y), by using the estimate 
¥(A f]B) < y/F(A)F(B), we have 

(it) < I EWR(B Sl - B S3 )y/R(B S2 - B Si )}ds, (3.9) 

J St 

(Hi) < [ E{ v / R(B Sl - B Si )^/R(B S2 - B s . A )}ds. (3.10) 

J s t 

Since R(x) = Yli=i(^~ \ x i\)^\x t \<i is bounded and compactly supported, we can find / £ 5(M rf ) such 
that y/R can be replaced by / in the above equalities. Writing by Fourier transform f(B Si — B s . ) = 

(2^ k*e* {B ' t ~ B ' i) f(Z)<%, we have 

^(U) < J^- d [ [ e-|l« 1 l 2u2 'e-5l« 1 +^l 2 -3* e -|l6l 2 ^/ (6) / fe)deide2dU) (311) 

i(«0 ^ T^W / / e-4"^l a '**e-il« 1+ &l au »'e-J^l a «»'/«i)/(6)deidfedtt, (3-12) 

¥ (2ir) M J R2d J Ui >o,Y? i=1 ui<i 



which leads to 

5S(»0 ^ 7A2S [ [ e-^+M 2 "*' I p K&K&dtid&duidua, (3.13) 

S(*») ^ ?Ad/ / e -|iaiVt 1 /(^ 1 )/( fe )de 1 dg 2 du 1 du 4 .(3.14) 

^ (2vr)^ 7 R2d 7 Ui >o, Ul +« 4 <i Ifillfi + 61 

By Lebesgue dominated convergence theorem, we see that as t — )• oo, (ii)/t 2 — > and (iii)/t 2 — > 0. 
For (i), by change of variables and integrating in U, we have 

W = E/ / Qs2-si(y2-yi)qs3-s 2 (y3-y2 + k)q S4 - S3 (y4-y3)dyds, (3.15) 

* J St Jl4 o 

where Jo = ( — 5' 3) ■ We ^ rs ^ integrate 2/1,2/4, write them in the Fourier domain and denote 
e(£) = ^{l/,,} = rifc=i S ™ fc 1 then change variables u; = Si — Sj_i,i = 1,2,4, U3 = (53 — 82)/*, so 

fc ^ ' ■/tti>0,U3<l,Ul+tt2-Hl4<t(l-U3) J/g ^K 4d 

<7t« 3 (y3 - 2/2 + k)dy 2 dy- i d^idi A du. 

(3.16) 
Integrating in m, U2, M4, we have 

- / e^5l?il 2 «2 e "2l€4l 2 «4 rfnidu2dn4 



ab \ Tb(b — a) Ta(a — b) 

ft- i 

2IS1I )° — 2 



if let u = u 3 ,a = I|^| 2 ,ft = ±\U\ 2 ,T = t(l - u 



We point out that rb( .^_ n (1 — e bT ) + Ta ,_ h - ) (1 — e aT ) is uniformly bounded by Lemma 

5.1. Therefore, \ f u ^ Ul+U2+Ui < t{l _ u) e^^eT^I 2 -^^^ = |^(1 + (*)) with (*) 
uniformly bounded and going to zero as t — > 00. Now we have 

@ = E7^2l/ / / ^2l^ ( 1 +(*)) a (^) g (^)e" ift3/2 e~^ 3 g to (2/3-2/2+fc)d2/2^ 3 ^i^4^. 

(3.18) 
Now J R2d J2 k Qtu(y3 - 2/2 + fe)l/ (2/2)1/0 (2/3 )e 4§lS/2 e ^ 4V3 dy 2 dy 3 is uniformly bounded since 



X M E (2 7 rfa)^ e " ly2 " y3 " fc|2/2falfo(y2)ljo(y3)e "^ 12/2e "^ 2/3dy2 ^ 31 



jR2d ^ [2iTtu) d i z ^7[-i,i]d [2-Ktu) a l z 

(3.19) 
In addition, , f rjrg re e(£i)e(£4) is integrable, so we have 

n) r(l+(*))e(6)e(^)e-^ y2 e-^ 3 ^g in (2/3-2/2+fc)ci2/2d2/3^i^4dn. 



i 2 (2vr) 2rf y u6{0il) y 7o2 y R 4 d |6i 2 iC4| 

(3.20) 



On the other hand, for fixed u £ (0, 1), we have as t — > oo 

f R2d E (^e-l,3-. + ^/- 1/o(y3)1/o( , 2)e -^, 2e -^^ y2dy3 

k y ' (3.21) 

-> / l/ (y 2 )l/o(y3)e-^ 1J ' 2 e-^^dy 2 dy3 = £(&)£(&), 

jR 2d 

and thus by Lebesgue dominated convergence theorem, 
To sum up, 



2 ^ 2 



E{£>(*) 2 ) a } = 8 (« + («) + («0) ~ ^55 (£ f^) = ^ (3 " 23) 

as t — )• oo. The proof is complete. □ 

Remark 3.4. We mention that if V(x) is a mean-zero stationary Gaussian field with fast-decaying 
covariance function R G 5(IR d ), and X(£) = J* V(.B s )ds, then E{X(t) 4 } < Ct 2 for some uniform 
constant C. We will need it in the proof of tightness. To see this, we write 

^E{X(t) 4 } = § / V E{R(B n - B T2 )R{B T , - B T4 )}ds, (3.24) 

* * J0< Sl <s 2 <s s <s 4 <t {n}={si} 

and the terms containing R(B Sl — B S3 )R(B S2 — B S4 ) or R(B Sl — B S4 )R(B S2 — B S3 ) go to zero as 
t — > oo by the same proof as in Lemma 3.3. For the term containing R(B Sl — B S2 )R{B SZ — B Si ), 
we do the same calculation, and get 

(0 = 72TTx2d I I e-^l 2 ^-^e-^l 2 ^-^i?(6)%2)^6^ (3.25) 

First integrate in s, then let t — > oo, we will see that (i) — >• 2 f , 2 \ d J Rd "ra^^C) • Therefore, we 
have K{X(t) 4 } < Ct 2 for some uniform constant C. 

From Lemma 3.2 and 3.3, we know that YlkPt{k) 2 /t — > a 2 d in L 2 , i.e., the variance of X £ {t) 
conditioning on B, U converges in 1? . 

Now we can prove the convergence of finite dimensional distributions, i.e., V0 = to < t\ < ... < 
£jv < T, aj E E fixed, we have 

N N 

E{exp(t^a i (X e (t i ) - X £ (ij-i)))} -> exp(--a^a|(£ J - i,-_i)) (3.26) 



2 
i=i i=i 



as e — )• 0. 



Actually, we have 

N N 

^2otj(X e (tj) - X e (tj-x)) =J2a j eJ2^ k (p t:j/£ 2(k) - p tjl/£ 2{k)) 

3=1 N k (3-27) 

= ^2^^2aje(p tj/£ 2(k) - p tjl/£ 2(k)). 
k j=i 

For sum of independent random variables, central limit theorem comes from the following two 
lemmas by the Lindeberg condition [5]. 

Lemma 3.5. J2k \E,j=l a j £ (Pt 3 /e 2 ( k ) ~ Ptj^/e 2 ^))) ~> °~d EjLi "jfe ~ *i-i) in probability. 
Proof. We write Y^k [JLjLi a j £ (Pt J /e 2 ( k ) ~ Pt,_i/ E 2 (*0)) = (*) + (*«)» where 

N 

(*) = ^«^ 2 ^(^./ e2 (A:)-p tj _ l/£2 (/c)) 2 , (3.28) 

j=i ft 

(**) = X] «i«i £2 X)^*iA 2 ( fc ) _ ^-iA 2 ( fe ))(P%A 2 ( fc ) _ P%-i/£ 2 ( fe ))- ( 3 - 29 ) 

i,j=l,...N;i^j k 

By stationarity and the fact that (^2 k Pt(k) 2 ) 2 ft — > a\ in L 2 we know (i) — > a 2 , Ylj=i a2 (*j ~~ tj-i) 
in 1? . So we only have to show (ii) — > in probability. 
For any i > j, consider 

nY^iptM k )-pn-t/e<k)){p tJ i £ 2{k)-p t ._ 1/E ,{k))}=nr ' v(B s )ds V ' v(B s )ds}, 

k Jti-i/e 1 J*i-i/e 2 

(3.30) 

by similar calculation as in the proof of Lemma 3.2, if denote S± = tj — tj—1,82 = U-\ — tj,8$ = 
ti — ti-i, we have 



U/s 2 rtj/e 2 

s 2 E{ I V{B s )ds / V(B s )ds} 



k-i/e 



(2vr) d J Rd |£| 



.2 



£ I _1|C|2*2 l|g|2 J1+J2 1 |e|2 S2+S3 1 |g|2 Jl±J2±J3 , 



(3.31) 



and since <5j > 0, by Lebesgue dominated convergence theorem, we have E{e X^fcO-We 2 ^) ~~ 
PU- 1 /e 2 (k))(Pt/e 2 (k) ~ Pt _i/e 2 (^))} "" ^ 0; hence (m) — >• in probability. The proof is complete. 
D 

Lemma 3.6. sup fc e^/ e 2 (k) — > m probability. 



10 



Proof. We calculate the n— th moment, E{| sup k p t / £ 2(k)\ n } < ^2 k E{p t / £ 2(k) n }. By similar calcu- 
lation as in the proof of Lemma 3.3, we have 



y2E{p t/e 2(k) n } = n! / / T qsh-sk-Axk ~ x k _i)dxds 

I J0<s 1 <s 2 <...<s n <t/e 2 Jl% k ^ n 

t ( f^ e2 f 

— n '~2 / / Qu(x)dxdu 



' '?„-, (3.32) 



with Iq = (—5, \) d , and we check that f 2I q u {x)dx < lo<«<i + , 2 \d/2 ^ d ^-u>i- Because d > 3, 

we have e 2 J2k^{Pt/e 2 (^) n } ^ s uniformly bounded. So P(sup fc ep t / e 2 (fe) > 5) < ^e n ~ 2 for some 
constant C. Pick n > 3, the proof is complete. □ 

3.1.2 Tightness 



Proposition 3.7. X £ (t) = e f^ 2 V{B s )ds is tight in C([0,T]). 

Proof. Consider X £ (t) - X £ (s) = e fj,^ V{B s )ds. We claim that E{(X e (t) - X £ (s)) 4 } < C(t - s) 2 
for some uniform constant C. Then tightness follows from [4] given the fact X E (0) = and 

F(\X £ (t)-X £ (s)\>X)<^(t-s) 2 . (3.33) 

To prove the claim, we need to bound the fourth moment of V(x). By Lemma 5.2, we have 

E{(X £ (t) - X £ (s))*} <Ce 4 f Yl E {^ 2 (\B T1 - B T2 \)^(\B T3 - B T4 \)}ds 

J0<si<S2<S3<S4,<(t-s)/e 2 r T .i_r g .i 

(3.34) 
for some bounded and compactly supported function eft. By Remark 3.4, the proof is complete. □ 



3.2 Convergence to a homogenized equation 

Before proving the main theorem, we point out that without much extra effort, Theorem 3.1 can 
be extended as follows. 

Proposition 3.8. Define X £ (t) = | J Q (V(-f-) — V(-f-)) ds, where Bj,B 2 are two independent 
Brownian motions, then 

X £ (t) =► V2a d Wt (3.35) 

weakly in C([0,T]), where Wt is a standard Brownian motion. 

It is straightforward to check that the proof of tightness does not change since we observe that 

n\X £ {t)-X £ {s)\ > \)<F(\X l£ (t)-X l£ (s)\ > ±)+F(\X 2e (t)-X 2e (a)\ >^)< ^(t-s) 2 , (3.36) 
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HX i£ (t):=y*V(^)ds. 

For convergence of finite dimensional distributions, since X £ (t) has the same distribution as 
£ Ylk^k{p\/ e 2{k) — p 2 / £2 {k)), where p\{k) = J* l[.B»+m=fccfe, we need to prove Lemma 3.5 and 3.6 

with p t / £ 2 (k) replaced by p\/ e 2 (k) — p\i £ i (&) . The proof of sup fc £\p]/ £ 2 (k) — p\i e i \ — > in probability- 
is trivial given Lemma 3.6, so we only need to check the proof of Lemma 3.5. 

Let X £ (t) = e Jq /£ \v(BI) - V{B 2 s ))ds. We look at J2f=i a j{ x Sj)- X e{tj-i)), where = t < 
ti < ... < tjv < T, which could be written as 

N I N N 

Y,<Xj{ X e(ti)-Xe(ti-l)) = E&M E a >t/^ -p\-,/e^ k )) " X>i(P? i/e a(*) ~ ^_ l/e *(*)) 
3=1 k \j=l 3=1 

(3.37) 
Our goal is to show that £ fe e 2 (EJLi a j (pj i/e a (*) - Pt^i/e 2 ( fc )) ~ ££=l «j (Pt, /e a (A;) - ^,_ l/e a (*))) 

converges in probability to 2a 2 , X^ 7 =i °^j(tj ~ tj—i)- We only have to deal with the cross term here, 
i.e., 

N N 

\^Y,Y, a M ] /e^)-p\_ l/ Ak))Y, a ^l/e^ k )-Pl-M k ))\ 

(3.38) 

N N K ' 

<^EEN^/-»( fc )-^-i/«-( fc ))EN(p?i/«-( fc )-^-i/«-( fc ))' 

k j=l j=l 

and the following lemma would conclude the proof of Proposition 3.8. 

Lemma 3.9. Vi, j, we have £ 2 E{J2 k {pl /£2 {k) -pl_ i/£ 2(k))(p 2 j/£2 (k) - p 2 t] _ i/£2 (k))} -> as e -> 0. 

Proof. We have 

' £ 2 (e -i|^_ l/£ 2 _ e-ifclV^e-^I'Vi/* 3 _ e -§KI V* 2 ) J-%)d£ -> 



(3.39) 
as e — > by Lebesgue dominated convergence theorem. □ 

Remark 3.10. Note that the same result holds for X £ (t) := - L ( V(— s -) + V(— -) ) ds. The proof 
is the same. 

Now we are ready to prove our main theorem in the short-range correlation case. 

Theorem 3.11. In the short-range correlation case, we have ¥,{\u £ (t,x) — Uo(t, x)\ 2 } — > as 
e — > 0, where u £ and uq solve the following PDE, respectively, given the same initial condition 
u £ (0,x) = u (0,x) = f{x) for f E C b (R d ): 

d t u £ = -Au £ + i-V(-)u £ , (3.40) 

d t u = -Au - -o- 2 d u Q . (3.41) 
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Remark 3.12. If we further assume / E L 1 (M rf ), since the solutions are both bounded by <!>, which 
is in L 2 (R d ) and solves ® t = ^A3> with initial condition <&(0,x) = |/(x)|, then J Rd E{\u £ (t,x) - 
uo(t, x)\ 2 }dx — > as e — > 0. 

Proof. By the Feynman-Kac formula, we write u e (t,x) = E B {f(x + B t )exp(ij j Q V ( x+ £ Bs )ds)} . 
We claim that {B t , \ Jj V(^^)ds) => (iVi, iV 2 ) in distribution as e -> with JVi ~ iV(0, £), iV 2 ~ 
N(0, a\t) and iVi, iV 2 being independent. To see this, since V is stationary, we only have to look at 

{eB t / e 2 , e f V(B s )ds) by scaling, and if we compute the characteristic function by conditioning 
on B, U, we have 

~t/e 2 pt/e 2 



pi/e- ft/e- 

E{exp(i/9ie5 t/£ 2 + 9 2 e / V(B s )ds)} =E{exp(i9 1 £B t/s 2)E{exp(ie 2 £ / V(B s )ds)\B,U}} 
Jo Jo 

ft/e 2 i i 

E{exp(i0 ie 5 t/£2 )E{exp(i0 2 £ / V(B s )ds) - exp{- -a 2 d e^t)\B , U}} + E{exp(i9 1 £B t/£ 2)exp(--a 2 d 9 2 2 t)} 

Jo l l 



'0 

(3.42) 
and (ii) = exp(— i# 2 £ - \9\a\t). Consider 



E{exp(i0 2 e I ' V(B s )ds)-exp(-^a 2 d 9 2 t)\B,U} = E{exp(i9 2 Y,^kPt/ £ 2(k))-exp(-^a 2 9 2 t)\B,U}. 

k 

(3.43) 
By Lemma 3.3 and 3.6, for any sequence e, — > 0, we can find a subsequence, still denoted by £j, 
such that zlYlkPt/e 2 ^) 2 — > o\t and sup fc £j£> t / £ ? (k) — > almost surely. So by the central limit 
theorem, Ylik £ ^kPt/£ 2 {k) =^ iV(0, <r^t) when freezing B,U. Thus we see (i) — > 0, and the claim is 
proved. 

Now we have 

rt/e 2 

E{u £ (t, x)} = E{f(x + eB t/£2 ) exp(i£ / V(B s )ds)} 

Jo 

^E{f(x + JVi) exp(ziV 2 )} = E{/(x + iYOlEjexp^iV,)} = E B {/(x + B t )} exp(-^t) = u (t, x). 



(3.44) 



;1 rVT/^ + B ^ TA/ ,:r + - B ?' 



Next we consider \u £ (t,x)\ 2 = E B {f(x + B})f(x + B 2 )exp{i £ -£(V(^f±)-V{^f^))ds)} for 
independent Brownian motions B ,B 2 . By the same discussion as before, we can show 

(B},B 2 ,- f\v(^^)-V(^^))ds)^(N 1 ,N 2 ,N 3 ) (3.45) 

£ J £ £ 

weakly as e — > 0, where iVi ~ N(0, t),N 2 ~ iV(0, t), N 3 ~ iV(0, 2a d t) and they are independent, so 

E{\u £ (t, x)\ 2 } = E{f(x + Bl)f(x + B 2 )exp{i- f (V{^±^-) - V(^±^))ds)} 
->E{f(x + N 1 )f(x + N 2 )exp(iN 3 )} = E{f(x + N 1 )}E{J(x~+N 2 )}E{exp(iN 3 )} (3 - 46) 



--E B {f(x + Bl)}E B {f{x + 5 2 )} exp(-^t) = \u (t, x) 
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so the proof is complete. □ 

In the short-range correlation setting, we have assumed that {£&} are i.i.d random variables, 
which simplifies the proof of central limit theorem. However, we expect certain mixing condition 
would suffice, i.e., Theorem 3.11 should hold when we only assume (2.9) for <f> satisfying some 
integrability condition. Note that the order by which different £& join the summation depends 
on the Brownian path, and when d > 3, Brownian motion is transient so statistically the mixing 
condition still holds for the summation of £& , but we need to quantify the transiency of Brownian 
motion, which is a non-trivial task not considered here. 

4 Long-range correlation 

We recall that in Theorem 3.11, the constant potential in the limiting equation is given by a^ = 
( ^3-r(| — 1) J Rd , iaJ 2 dx) so that as long as R(x)/\x\ d ~ 2 is integrable, a^ is well-defined. There- 
fore in the long-range correlation case 1, since R(x) ~ Cd|x|~ a with a > 2, we still expect a 
homogenization result to hold. 

In the following, we will in turn deal with the cases a E (2, d) and a E (0,2), where the sizes 
of the potential are chosen to be different, i.e., - and -^. We observe a transition occurring at 
a = 2. 

4.1 Case 1: a G (2, d) 

In this setting, the proof is similar to the short-range correlation case. Since {£&} are assumed to be 
Gaussian, we do not need to apply the central limit theorem and only have to show the convergence 
of the variance. We repeat what has been done in Lemma 3.2 and 3.3 with sometimes significantly 
more involved calculations. 

We write X e (t) = £ f V(B s )ds = £j2k£kPt/e 2 (k) with the same notation, with a variance 
conditioning on B, U given by 

E{X £ (tf\B,U} = e 2 ^r(m - n) Pt/e 2(m) Pt/e 2(n), (4.1) 

and the convergence in probability is given by the following two lemmas. 
Lemma 4.1. E-Em,n r ( m ~ n)p t {m)p t (n)} ~ a^t. 
Lemma 4.2. E{(^ mn r(m — n)p t (m)pt(n)) 2 } ~ a^t 2 . 

Proof of Lemma 4.1. Let X(t) = J V(B s )ds = J^k €kPt(k) so that E{^ mjn r(m-n)p t (m)p t (n)} = 
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E{X(t) 2 }. We have 

E{X(t) 2 } = I I E{R(B S - B u )}dsdu = 2 I I I R(x) ( J—Y exp(-^-)dxduds 
Jo Jo Jo Jo JR d \2ttuJ 2u 

1 1^ r I ^\^- 2 exp(-X)l, xl2 l 0<s<1 dXd S dx (4-2) 

Jo Jo Jm d \ x \ A> ^r 



vr d/2 Jo Jo Jw \ x \' 

«(*);__ t w<* ^ f R ( x ) 



-L f 00 At" 2 exp(-A)dA / ^dx = * r(£ - 1) / 

7T d/2 JO JVL*\A 7T ' 2 JM 



, \rl— 2 



<ix 



as t — > oo by Lebesgue dominated convergence theorem, since R{x)/\x\ d 2 is integrable. The proof 
is complete. □ 

Proof of Lemma J^.2. Since ^ are Gaussian, we have E{(X^m „ r(m—n)pt(m)pt(n)) 2 } = gE{X(t) 4 } 
for X(t) = J Q t V(B s )ds = EktkPt(k). Given that 

4 

nil V ^ = E r(h)r(k 2 )¥([ Xl + U}- [x 2 + U]=h, [x 3 + U\- [x A + U} = k 2 ) 

i=l ki,k2 

+ J2 r(k 1 )r(k 2 M[x 1 + U]- [x 3 + U} = h, [x 2 + U] - [x 4 + U] = k 2 ) (4.3) 

k\,k2 

+ E r(fei)r(fc 2 )P([xi + U\- [x 4 + U} = k u [x 2 + U]- [x 3 + U] = k 2 ), 

ki,k2 

we have E{X(t) 4 } = 24((i) + (ii) + (Hi)), where 

(i) = E r(k x )r(k 2 ) f V([B 31 +U}- [B S2 + U} = h, [B S3 + U\- [B Si +U] = k 2 )ds, (4.4) 

(ii) = E r(h)r(k 2 ) f F([B S1 + U] - [B s , +U} = h, [B S2 + U] - [B S4 +U] = k 2 )ds, (4.5) 
ki,k 2 Jst 

(Hi) = Y r(h)r(k 2 ) i F([B S1 + U] - [B Si +U] = k h [B S2 + U] - [B S3 +U] = k 2 )ds (4.6) 
fei,fe 2 St 

for St = {0 < si < s 2 < S3 < 54 < £}. Next, we deal with them in turn. 



For (i), after some calculation and change of variables, and let Iq = (— A, A) , we have 



CO j. 

t 2 t 2 



n.i,k 2 ,k^ St ! o 

4^ 7,*, L 1/o(xi)lj ° (xi + y)^°(^)iio(^ - *) E ^ ^d-2 E 



fel &2 



d rt d 



Ew( x 4 -2- xi -y + fc)A| e A2 A 4 2 e A H , [ M+fcl |2 lz+k2l 2 ^^dydzdx 1 dx 4 du 1 du 3 dX 2 dX 4: . 

J—* |Ml+«3+ 2Aat + 2A4t S-Lj 

(4.7) 
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We note that 

|l, (x 1+y )l /() (x 4 -.) ^ +fc|J _ 2 L k + fc2 |^- 2 l < E W^) , +fe|d - 2 E lM °( g ) jI+fap^ ' 

^£7 ' J_ i ■ £1 \ i t7-L I ' it \ 

1 K 2 fel K 2 

(4-8) 
and since / Rd ^fc 1 2J (^) ^JpffiU ^ = / K d j^2 Efc 1 2/ (2/— fc ) k(^) 1^2/, by the fact that \r(k)\ < 
\k\~ a , we have ^2 k 1 2/0(2/ ~~ &)|r(A;)|dy < 1 A |y| _a , thus using the fact that a > 2 we know 
Sfe l2Jo(2/) 1 +Ld~2 i s integrable. On the other hand, for fixed u% € (0, 1), we have J2 k q U3 t(%4: — 
z — xx — y+k) — > 1 as t ->■ 00. Moreover, J" R2d l/ (xi)l/ (x 4 ) ^ fe q^ix^- z - x\- y + k)dxidxi 
is uniformly bounded. In summary, we can apply Lebesgue dominated convergence theorem 
to get 



^ L L l7 « (Xl)1/0(Xl + V)l*(^)l/o(^ - *) E ly + frV-2 E ' 



r 47r JR« JR*_ ^iy + K il ~J~ 1^ + ^21 



A| e~ 2 A| e~ 4 l{ Ul+U3<1 }dy<izdxidx4ci'Ui<iM3GL\2(iA4. 



We check that 



ms^^-d'O^^ww+ioE 



r(fc) 

W+k\ 

— d n d --i) 2 ([ ^ 



^ d _ 2 dxdy 

2 



(4.9) 



(4.10) 



(ii). By a similar calculation, we have 

l ( ' ~\ l 1 r r 

— — =|— E / / r(/ci)r(fe 2 )fe-si(a;2 - »i + /c4-/ci)g s S2 (x 3 - x 2 - fc^-ssO^ - x 3 + /e + k 2 )dxds\ 
1 * kfto,k Js * Jl & 

^ v^ f f \ r (ki — k)\ Ir(/c2 — fc) I 

~ Z^ / / 1 2/o( a; ) 1 2/o(2/) 1 2/o(^)| — 1 , | rf _2 i — 1 1. irf- 2 g fa 3( a: ~ k)l ul+U3 <idxdydzduidu 3 . 

kl .k 2 ,k jR3djR l \y + ki\ \z + k 2 \ 

(4.11) 
Define G{y-k) = J2 n hi a (y~ k-n)\r(n)\ and F(k) = j Rd r^=^G{y -k)dy. It is straightfor- 
ward to check that G(y—k) < lA|y— &| _Q , and thus by Lemma 5.3, we have F(k) < lA|/c| _a+2 , 
which leads to 



— TT~ ~E / / 1 2/o( a: )( 1 A |Ma-2 ) 2 g*"3( :c ~ fc)l Ml+M3 <i(ixdtti(JU3 

= / / qtu 3 (x)Y]l2i (x - k)(l A —^ ZI )dxl Ul+U3<1 du 1 du 3 . 

Jul Jm. d u \ k \ 



Since Y^k^i { x — k)(l A pU) < 1 A \x\ 2a + 4 j by Lemma 5.4 and Lebesgue dominated 
convergence theorem, '-^- — > as t — > 00. 
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{in). Similarly, we have 
{in) 1 



+2 X) / r ( k i)r(k 2 )q S2 - Sl (x 2 -xi + k + ki)q S3 - S2 (x 3 - x 2 -k 2 )q S4 - S3 (x A -x 3 -k + k 2 )dxds 
ki,k 2 ,k o 

+2 X / r ( fcl ~~ k2 ~ k ) r ( k 2)q S2 - Sl (x2 -xi + ki)q S3 _ S2 (x 3 -x 2 - k 2 )q S/t - S3 (x 4: - x 3 - k)dxds, 



ki,k2,k o 

(4.13) 
and with a change of variables, 

|(m)| v-^ /" /" ,, , /nK*i-*2-*)I l r ( fc 2)| 



~ X / / 1 2/o( 2; ) 1 2/o(y) 1 2/o(2) — p-T^-i T ^-j Z Kq t u i {z-k)l Ul+u . i <idxdydzduidu 3 . 

,. ,. i. JR 3d JR 2 , \X -\- K\\ \y — «2| 

(4.14) 



i 2 

ki,k2,k 



First of all, by the same discussion as for (ii), we have 

jL E w*) 1 ^^ 1 ^ ~ x A " 2 + "" Q - (415) 

fei 

Next, we need to estimate (*) = f Rd Y, k2 hi (y) fi^ffe? (1 A | fc2+ ^-0 ^> and we write (*) = 
/ R d ]^F(y,k)dy for 

F(y, fc) = £ hi (y ~ fe)(l A -^-)(1 A |fc _L_ 2 )- (4-16) 

Since F(y,k) < 1 A |y| _c \ (*) is uniformly bounded. Now we consider the case when \k\ is 
sufficiently large, and write (*) = (1) + (2) with 

J\v\<m \y\ k2 |« K 2\ 

(2) ~ [^ r \^^ l2Io{y ~ k2) V^ {1A \k-l\^ )dy (418) 

J\y\>M \y\ k2 \ K 2\ \K K 2 \ 

for some fixed constant M. Note that for (1) we have \k 2 \ < M + \fd, so (1) < it-ia-n ■ F° r (2) ; 
we have £ fca l 2Io (y - k 2 )^{\ A |fc _ fc ^ |Q - 2 ) < J2k 2 W^ ~ fc 2)^(l A \ k _ k \\ a - 2 ), and thus 

(2) ~ 4m m^ £ Wy " * " " 2)(1 A w^- (419) 

Let G(y - k) = £ fc2 l 2/o (y - k - k 2 )(l A p^U) < 1 A \y - k\ 2 ~ a . Therefore, 

( 2 )^/ r^W 1A i^ \ a ~2 )dy- ( 4 - 2 °) 

J\y\>M \y\ d I+a \k-y\ a 2 

Note that the RHS does not fall into the form of integral in Lemma 5.3, but following the 
same prove as in the lemma, we can show (2) < ,, , a _ 2 . To sum up, we have shown that 

(*) < 1 A jj^i, so J Rd J R 2 Ylk l 2i {z){*)qtu z {z - k)l ur+U3 <idzdu\duz -> as t ->■ oo. 
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The proof is complete. □ 

From Lemma 4.1 and 4.2, we have 

E{X £ (t) 2 \B,U} = e 2 ^r(m - n) Pt/e 2(m) Pt/e 2(n) ^ a 2 d t (4.21) 

m,n 

in L 2 as e — > 0. We claim that X £ (t) converges weakly to some Gaussian iV(0, o~ d t). To see this, we 
only have to note that 

E | e iex £ (t)| = E/]E/g» ee £fe£kP*/ e a( fe )|£ ) [m = E{ e 4 j2£2 L/( m -")p«/ E 2( m )! , i/ E 2(")} i (4.22) 

Before proving the main theorem, we show the following lemma, which basically says that the 
interaction between two independent Brownian motions in a random scenery is small. 

Lemma 4.3. For any function R(x) such that R(x)/\x\ d ~ 2 is integrable, we have 

-f I R{B l (s)-B 2 (u))dsdu^0 (4.23) 

t Jo Jo 

in probability as t — > oo, where B l ,B 2 are two independent Brownian motions. 
Proof. In the proof we could assume R{x) > and by direct calculation, we have 

-/ E{R(B 1 (s)-B 2 {u))}dsdu = - / / R{x)q s+U {x)dxdsdu 
t Jo Jo * Jo Jo Jm. d 

= 7 / / — Jni^2 A ^~ 2ex P(~ A ) 1 (0,*)( u ) 1 r Ml. \,\ 2 S X ) dXdxdu (4.24) 

T JR\ JR d 27T2 \ x \ ( 2(t+«) ' 2u I 

1 R(x) d 

r A" 2 exp(-A)l (oa) (ti)l N 2 M 2 (X)d\dxdu, 



'/?2_ JR d 27T2 Fl "> 2t(l+u) ' 2tu > 

and now we only need the Lebesgue dominated convergence theorem to finish the proof. □ 

Now we can prove the main theorem for this case. 

Theorem 4.4. In the long-range correlation case 1, we have K{\u £ (t,x) — Uo(t, x)\ 2 } — > as 
e — > 0, where u £ , uq solves the following PDE respectively given the same initial condition u E (0, x) = 
uo(0,x) = f(x)forfeC b (M d ): 

d t u e = \&u £ + i-V{-)u e , (4.25) 

d t u = -Am - -<y 2 d UQ. (4.26) 



Remark 4.5. Similar to the short-range correlation case, we have J Kd E{|it e (£, x) — Uo(t, x)\ 2 }dx — > 
ase^Oif / € L x (R d ). 
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Proof. The proof is similar with what has been done for Theorem 3.11, so we do not present all 
the details here. 

First of all, we show the convergence of K{u £ (t,x)}. Since X e (t) = e J V(B s )ds weakly 
converges to some Gaussian N(0,a^t), we only have to show that eB t / £ 2 and X £ (t) are asymptot- 
ically independent, and this is the same as in Theorem 3.11. Therefore, we have E{u £ (t,x)} — > 
E B {f(x + B t )} exp{-±a 2 t) = u (t,x). 

Secondly we prove the convergence of 

rt/e 2 

K{\u £ (t,x)\ 2 } =E{f(x + eBl /£2 )f(x + eB^ /£2 )eMieJ (V(Bl) -V(B 2 ))ds)} (4.27) 

for independent Brownian motions B l , B 2 . We point out that e L {V{B\ ) — V(B 2 ))ds converges 
weakly to some Gaussian N(0, 2o~ 2 j t), which comes from the fact that £& are Gaussian and Lemma 
4.1, 4.2, 4.3. The rest of the proof is the same as in Theorem 3.11. 

The proof is complete. □ 



4.2 Case 2: a e (0,2) 

In this setting, we first present some properties of the random field. Since $ satisfies E{<&(r/fc) 2 } < 
oo, we can write £& = &(r]k) in Hermite polynomial expansion 

oo 

& = *(Tfc) = Z)-Tflnfak), (4-28) 

n=0 

where H n (x) = (— l) n exp(^-)^: exp(— ^-) is the n— th order Hermite polynomial and V n = 
K{H n (rjk)^(r]k)} ■ By the assumption that Vq = 0, Vi / 0, we have 

oo y 2 

n=0 H ' (4.29) 

Since Yl^=o ~nf < °°j we see that r(k) ~ V 2 p(k) as \k\ — > oo. By (2.6) and the fact that p{k) ~ 
c <i ]T[i=i l^i| -a S we know that R(x) ~ V^c^ ]^[ i=1 |cci| — Qi as min i=1) rf|xj| — > oo. Because a = 
Yli=i a i < 2, the correlation function R(x) decays sufficiently slowly so that the integral X £ (t) = 
—^j2 J V(- £ ± )ds is able to memorize the Brownian path in the limit e — > 0. 

The assumption that V\ ^ is crucial for the appearance of Gaussian white noise in the limiting 
equation, and it turns out that we can reduce it to the case $>(x) = x, namely V(x) = rj\ [X+u ^ so 
conditioning on B, U, X £ (t) = -h$ f V{-f)ds is Gaussian, and we can prove its weak convergence 
by calculating the variance. Before doing that, following [10] we define the solution to the limiting 
SPDE. 
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4.2.1 limiting SPDE 

We first define the formally-written random variable L J RtJ S(x — B s )W{dx)ds, where W(dx) is 
the generalized Gaussian random field independent from Brownian motion Bt, and the covariance 
function E{W (dx)W (dy)} = =a — r ^r^xdy. 

lli = l \ x i~yi\ % 

Proposition 4.6. Assume Yli=i a i < % and define Y £ (t) = J Q j Rd q £ (x — B s )W{dx)ds. Then Y £ (t) 
converges in L 2 as e — > to some random variable Y{t), denoted as Y(t) = L J Rd 5(x—B s )W(dx)ds. 
When conditioning on B, then Yt is a Gaussian random variable with zero mean and variance 



E{Y(t) 2 \B} = I I — — -i dsdu, (4.30) 

Jo Jo Ht =1 \Bi(s) - Bi(u)\<*< 

where -Bj(s) denotes the i—th component of B s . 

Proof. We first point out that the RHS of (4.30) is almost surely finite, and this comes from the 
fact that Yli=i a « < 2 and 

E{ I I —, dsdu} = I f —. dsduf]E{\N\- a '}, (4.31) 

Vo Jo nLi\Bi(s)-Bi(u)\°* Jo Jo | a -u|Eti«*/2 f = \ 

where N ~ N(0, 1). Secondly, we calculate 

^{Ye\t)}= I I I E{q s (x - B s )q E (y - B u )} * —dxdydsdu. (4.32) 

Jo Jo Jr™ IIt=i ki-wh 

By Lemma 5.5, J R2d q e (x - B s )q e (y - B u ) ^-^—^- dxdy -)• ^- ^.^^ as e -> 0. By 

Lemma 5.6, i.e., the inequality J R2 q e ( x i + Vi)Qe( x 2 + 2/2) 1 2/1 - y2\~ a dy±dy 2 < C\x\ - x 2 \~ a as 
long as a E (0,1), and by the Lebesgue dominated convergence theorem, we have the conver- 
gence E{Y 2 (t)} -> ft ^{ ^-j^^-^-^ jdsdu. Similarly, we can show E{Y £l (t)Y £2 (t)} -> 

L LE{— 3 — - 1 — }dsdu as £i,£2 — ► 0. Thus, we have shown that {Y £ (i)} is a Cauchy 

lli=l \Bi(s)— Bi(u)\ l 

sequence in L 2 since lim ei)e2 ^oIE{(Y ei (£) — Y £2 (t)) 2 } = 0. The limit is then denoted as Y(t) = 
Jq J^d8{x ~ B s )W{dx)ds. 

Next, we consider the conditional distribution. Since Y £ (t) — > Y(t) in L 2 , there exists a subse- 
quence £fc such that Y £k (t) —> Y(t) almost surely. Note that W(dx) and B t are independent, so the 
probability space is the product space. Then we know that conditioning on the Brownian motion, 
Y £k (t) —7- Y(t) almost surely as k — > 00, and this leads to convergence in distribution. Given B, 
Y £ (t) is Gaussian with variance 

r-t rt [■ 1 

-dxdydsdu 



E{Y 2 (t)\B}=[ f f q e {x-B s )q £ (y-B u ) l 

Jo Jo Jm. 2d Y\i=i \ x i 



yir 



'0 J° I\i=l \ B i( s ) - B i( u 

The proof is complete. □ 
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(4.33) 
dsdu. 



Remark 4.7. If we define Y l (t) = j Q j Rd 5(x — B l s )W{dx)ds for independent Brownian motions 
-B 1 ,!? 2 , the same proof implies that y 1 (t),y 2 (t) are jointly Gaussian with covariance function 
given by E{Y l (t)Y 2 (t)\B} = £ £ ^-^___ dsdu when conditioning on B\B 2 . 

Remark 4.8. Prom the proof of Proposition 4.6, we see that the distribution of J j Rd 5(x — 
B s )W{dx)ds does not depend on the starting point of Brownian motion. 



With random variable J J Kd 5(x — B s )W(dx)ds, we can formally write the solution to the SPDE 

d t u = -Au + iWu (4.34) 

with initial condition u(0, x) = f(x) by Feynman-Kac formula as 

u(t,x)=E B {f(x + B t )exp(i I [ 5{y-x-B s )W{dy)ds)}. (4.35) 

JO JR d 

We point that the u(t, x) defined as above coincides with the usual definition of weak solution to 
SPDE (4.34): 

Definition 4.9. A random field u(t,x) is a weak solution to (4.34) if for any C°° function cj> with 
compact support we have 

/ u(t,x)cf)(x)dx = / f(x)cj)(x)dx + - / / u(s,x)A(j)(x)dxds + i / / u(s,x)<f)(x)W(dx)ds. 

JR d JlRd 2 Jq J R d J J R d 

(4.36) 
Proposition 4.10. // E i=i a i < 2, u(t,x) is the weak solution to (4.34). 

The proof is a direct adaption of Theorem 4.3 in [10], and we do not present it here. 

4.2.2 Convergence to a stochastic equation 

First of all we write X £ (t) = -^ f V(—r)ds = -^ Ylk ^{Vk)Pt{k), where we have defined Pt(k) = 
Jo ^-\S^ +U i =k ds = f lB g ee(i k ~u)ds and l k = xf =1 [h, ki + 1). We reduce it to the Gaussian case by 
the following lemma: 

Lemma 4.11. Y e (t) := -^ X^fc(^K%) ~~ ViVk)Pt(k) — > in probability as e — > 0. 
Proof. Conditioning on B,U, we have 

oo T/2 



1 V 2 

nYe(t) 2 \B, U}=-^2Y1 -jP( m " n) k Pt(m)pt(n) 



m,n k=2 

c 

-~^ZP^ m ~ n) 2 pt(m)p t (n) 



(4.37) 



for some constant C, and the RHS can be written as 



RHS = — f f p £ { Bs Bu ) 2 dsdu (4.38) 



£ Jo Jo 
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if we define p e (|(x — y)) = /5(7£(to — n)) if x G e(I m — U), y G e(/ n — U). Since p is bounded, and 
by Lemma 4.13, we have 



RHS <C sup \p{k)\ / / — — l\ Bs -B u \>Medsdu 

\k\>M-Vd J o Jo [[ i=1 \Bi(s) - Bi(u)\ a * 

C /"' /"* 

+ 7^ / / l \Bs~B u \<Medsdu, 

e Jo Jo 



which leads to 

K{Y £ (t) 2 } < C sup \p(k)\ + -^ I f I F(\B 8 - B u \ < Me)dsdu. (4.40) 

\k\>M-Vd £ Jo Jo 

By Lemma 5.7, first let e — > 0, then M — > oo, the proof is complete. □ 

Now we can prove the weak convergence of X £ (t) . 
Proposition 4.12. X £ (t) => Vi^fcd J * f Rd 5(x - B s )W(dx)ds as e ->■ 0. 

Proo/. We write X £ (t) = ^ ^(^(j/fc) ~ Vir)k)pt(k) + -^ Y,k ViVkPt(k), and by Lemma 4.11, we 
only need to show the weak convergence of -^ Y^k Vir] k p t (k). 

By conditioning on B,U, we calculate the characteristic function 

1 V 2 6 2 

E{exp(i0^ ^ ViVkPt(k))} = E{exp(— ^- ^ P( m ~ n)p t {m)p t {n))} , (4.41) 

k m,n 

and claim that i J2m,n P( m ~ n)p t (m)p t (n) ->■ c d /<J jj ^ j^jp^)]^ dsc ^ almost surely. 
Actually, we write 

1 . 1 1 /"*/"' 1 

— YV-e(m - n))p t {m)p t {n) = — p £ (-(B s - B u ))dsdu (4.42) 

£ m,n £ E J JO £ 

where we have also defined p £ (^(x — y)) = p(^e(m — n)) if x G e(/ m — C), y E e(/n — ^0- Thus by 
our assumptions on p(fe), we see that p £ {-(x — y)) ~ c^e a LJi=i l x « ~~ 2/il"" 1 for all x,y satisfying 
Vi,Xj ^ yj. Now we only need a bound to pass to the limit, and this is given by Lemma 4.13. 
Therefore, we have 



E{exp(i0— *- J^Viri k pt(k))} ->E{exp(--9 2 V?c d I I —, dsdu)} 

=E{exp(i9Viy/c~i I I 5(x - B s )W(dx)ds)} 

Jo JlBL d 



(4.43) 



as e — > 0. □ 

Til 



Lemma 4.13. \p £ {-(x — y))\ < Ce a rii=i \ x i ~ Di\ ai f or some uniform constant C independent of 
e and x, y. 
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Proof. Pick M > 0,x,y fixed. Let A = {i : \x{ — yi\ < Me}, since p £ {j(x — y)) = p(je(m — n)), if 
i g" A, i.e., \xi — yi\ > Me, we have \rrii — m\ > M — 1. So 

1 -^-r p a * -^-^ r Qi -^- It-— ii-\ ai 1 



k(-(a?-sO)| <cTT — -^ ^<cT~\ Tl —. — ^=cT\ 



A £ u *i\m i -n i \ lx i ^\\Xi-yi\-e\^ " ||£; - 2/i| - q Qi R ~ 2/in 

<c TT c a ' x 1 

" 11 (1 - 1/M)«i fa - yi\<* ' 

(4.44) 
On the other hand, for i £ A, we have riieA^I^* ~~ Ui\~ ai ^ FLeA M~ ai , so the proof is complete. 
□ 

Now we are ready to prove our main theorem. 

Theorem 4.14. In the long-range correlation case 2, we have u £ {t,x) =r- Uo(t,x) in distribution 
as e — > 0, where u £ ,uq solves the following PDE respectively given the same initial condition 
u £ (0,x) = uo(0,s) = f(x) for f G C b (M d ): 

d t u e = -Au £ + i— -fiV(-)u e , (4.45) 

2 e a i z e 

d t u = -Auq + iVx^rWuQ. (4.46) 



Proof. For fixed (t,x), we write 

Z £ :=u e {t,x) = E B {f(x + B t )exp(i^- f V(^^)ds)}, (4.47) 

Z :=n (t,x) = E B {f(x + B t )exp(iViyfe [ I 5{y - x - B s )W{dy)ds)}, (4.48) 

JO JR d 

and claim that Mm, n G N, E{Z £ m Z^} -> E{Z™Z*}}. Actually, we have 



J=l j=m+l j=l j=m+l 

(4.49) 
where i?j , J = 1, . . . , iV = m + n are independent Brownian motions. Since V(x) is stationary and 
all relevant functions are bounded and continuous, to prove the convergence of E{Z™Z™}, we only 
need to prove the convergence in distribution of 

N - 1 /•* R j 

W E :=J2^Bl+J2^^E y(-r)d* (4-50) 
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for dj,j3j £ R. We write W e = (i) + (ii) + (m) with 



AT 

J^aX, (4.51) 

i=i 

AT 



N i 

(m) = ^—^(^-^^(AO, ( 453 ) 

j=l £ k 

where we have defined pj-(k) = L 1 B j ds. (in) — > in probability by Lemma 4.11, and for 

[-^.+C7]=fc 

(i) + (ii), we calculate 

N 1*1 

E{exp(*0i(*) + *0 2 (u))} = E{exp(^ 1 ]Ta ? ^)e X p(--Vf0 2 2 ]T &/?,- £>(ro- n)rf(m)^(n))}, 

(4.54) 



0=1 i,j=l m,n 



and by the same proof as in Proposition 4.12, we have 



I ft ft 

— > p(m — n)pl(m)pl(n) — > / / — 3 : ds<iu (4.55) 



almost surely. Therefore, we see that 



N N 



(i) + (^) => V a^ + ^iv^ V Pj [ f 5(y-x- Bi)W(dy)ds (4.56) 

Pi Pi Jo J* d 

in distribution given Remark 4.7, so W e => Ylj=i a j^t +^iv / ^Si=i /^i Jo JiR d S(y — x — B J s )W(dy)ds 
in distribution. Thus the claim is proved. 

Note that \Z £ \,\Zq\ are uniformly bounded, if we let Z e = Z £j i + iZ £t 2,Zo = Zq^ + IZq^i 
the corresponding real and imaginary parts are uniformly bounded as well. From the fact that 
E{Z™Zf} -> EjZ^Z™}, we know \lm,n£ N, E{Z™ X Z£ 2 } -> EjZ^Z^}. So 

oo 1 

E{exp(^ 1 Z £il + i0 2 Z £>2 )} = ^ -E{(i^Z eil + i9 2 Z £)2 ) k } 

fc =° ; (4.57) 



J^ iE{(i0iZ o ,i + ^ 2 ^o, 2 ) fc } = E{exp(i0iZ o ,i + ^ 2 Z , 2 )}, 



fc! 

fc=0 



which completes the proof. 

□ 



24 



5 Conclusions and discussions 

In this paper, we have studied a type of partial differential equation with a large, highly oscillatory, 
random coefficient with either short- or long-range correlation in dimensions d > 3. Using a 
Feynman-Kac formula and weak convergence approach, we have shown that the solution to the 
PDE converged either to a homogenized or to a stochastic equation. The case d = 2 could be 
handled in the same way with a logarithmic correction showing up in the scaling coefficient (see, 
e.g., [7]). The reason for the dependence of the limiting equation on the correlation property of 
random potential could be explained heuristically as follows. The variance of the solution is affected 
by the interaction of two independent Brownian motions in the random scenery, i.e., L V{B\)ds 
and f V(B^)du, and we need to estimate the term f Q f Q R(B}. — B^)dsdu. If the random potential 
is short-range correlated, R(x) decays sufficiently fast, then the interaction is negligible in the 
asymptotic limit and the variance goes to zero. On the other hand, for the long-range correlated 
potential, this interaction remains non-trivial as e — > and the limiting equation remains stochastic. 

There are some possible directions to apply and extend our results. We mention three of them 
here. Firstly, we could consider the parabolic Anderson model ut = <^ n + ^e n ' m which case 
even weak solutions may no longer be well-defined. In the long-range correlation case, it could 
presumably be analyzed by the estimates in [10]. In the short-range correlation case, [2] proved a 
homogenization results for Gaussian potential when t £ (0, T) for sufficiently small T. It is not clear 
that for our potentials, the solution is well-defined for all t > 0. If we assume V(x) to be Gaussian 
as well, to recover part of the results in [2] , together with the weak convergence results given by [7] 
we still have to show the uniform integrability of exp(L V e {B s )ds) to pass to the limit. This could 
be done by writing the exponential function as a Talor series and estimate each term. We point 
out that in this case, diagrammatic expansions as done in [2] could be avoided. Secondly, in the 
long-range case, the condition V\ ^ ensures that the non-Gaussian case reduces to the Gaussian 
case so we have Gaussian noise in the limiting SPDE. If we still consider £& of the form ^(ijk) for 
some function $, or another type of random field V(x) = $>(v(x)) for some function <3? and Gaussian 
field v(x), the Hermite rank of $ determines the weak convergence limit of —^ J V(-£-)ds, and 
the case Vb = V± = 0, Vi 7^ would presumably produce some Rosenblatt-type noise in the limit, 
at least in the one-dimensional setting (see, e.g., [8]). Finally, in the intermediate settings, i.e., 
the long-range correlation case 1, we were only able to prove a homogenization result when £*. are 
Gaussian. At present, we cannot generalize our derivation to the case when £& = $(%), as what 
has been done in long-range correlation case 2. The reason is that the scaling 1/e makes all the 
terms in the Hermite expansion contribute in the limit and we cannot reduce it to the Gaussian 
case even assuming the Hermite rank of <I> to be 1. Another type of limit theorem for the sum of 
random variables with long-range correlation is needed. 
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Appendix: technical lemmas 

a (\ r - bT \ I b 



Lemma 5.1. T ui_ \ (1 — e bT ) + r a / a _j,\ (1 ~~ e aT ) *s uniformly bounded when a,b,T > 0. 



Proof. We only have to look at f(a, b) = ,,,"_ s (1 — e fc ) -\ — , _,-, (1 — e a ), and it is straightforward 

to check that for fixed a > 0, lim&_>.o f(a, b) = — 1, lim6_> a /(a, 6) = — — ~ 2+2e — is bounded, and 
lima^o /(a, b) = — 1. Assume there exists (a n ,b n ) such that /„ = f(a n ,b n ) — > oo, we know that 
there exists a subsequence a n — > oo, b n — > oo without changing the notation. If a n /b n is unbounded, 
we get the contradiction. So we assume that a n /b n — >• > 0, if ^ 1, we get the contradiction. So 
we assume that a n /b n — > 1, then we only have to write 

f, p — dn — p — On 

fn = (e— - 1)(- + i) + ^ e . _ 6 -► 0, (5.1) 

as n — )• oo and get the contradiction. Thus, the proof is complete. □ 

Lemma 5.2. Let V(x) be mean zero stationary random field with K{V(x) G } < oo, and satisfying 
strongly mixing property with mixing coefficient <f>, then we have 

\E{v( Xl )v(x 2 )V(x 3 )v(x 4 )}\<c J2 ^(|yi-y 2 |)^(|y3-y4|)Emx) 6 }l. (5.2) 

{yk}={x k } 

Proof. Let 2/1 and y 2 be two points in {xk}i<k<A such that d(yi,y2) > d(xi,Xj) for all 1 < 
i,j < 4 and such that d(yi, {2/3,2/4}) < d(y 2 ,{y 3 ,yi}), where {yfc}i<fc< 4 = {xk}i<k<i- We assume 
d(y3,yi) < d(jj4,yx). Therefore by (2.9), we have 

£ := Ie^si)^)^)^)}! < 0(2| yi -2/ 3 |)(iE{^(yi) 2 })^iE{(^(y 2 )^(y3)^(y4)) 2 })i (5.3) 

The last two terms are bounded by E{y(a;) 6 }6 and E{V(x) 6 }2 respectively. Because <p(r) is 
decaying in (0, oo), we have £ < (j)(\yi — 2/3|)E{V "(a;) 6 } 3 . On the other hand, if 7/4 is (one of) the 
closest point(s) to 2/2, the same argument shows that £ < 0(|j/2 — 2/4|)E-{V(x) 6 }3. Otherwise, 2/3 
is the closest point to y 2 , and we find £ < 0(2|j/2 — j/3|)1E{1 a (x) 6 } 3 . However, by construction, we 
have 

1 2/2 - 2/4 1 < I2/1 - 2/2I < I2/1 - 2/3 1 + 1 2/2 - 2/3I < 2 1 j/2 - 2/3 1, (5-4) 

so we still have £ < 4>{\y2 — 2/4|)E{V^(x) 6 }3 . To summarize, we have 

£ < <Pk\yi - Va\)<f>H\V2 ~ y 4 |)E{y(x) 6 }§, (5.5) 

and this completes the proof. □ 

Lemma 5.3. Assume ct,/3& (0, d) and a + /3 > d, then 

dy < C— Ij--, (5.6) 



R d 12/1 \ x — 2/ 1 



for constant C independent of x. Moreover, we have 



1 K 1 A rA^)^ ^ A wsW ( 5 - 7 ) 



Rd |2/|" F-2/1 
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Proof. Let p = \x\, and define B{z,r) = {x : \x — z\ < r}, (i) = {z : \z\ < \z — x\},(ii) = 
{z : \z\ > \z — x\}, we divide R d into three disjoint parts, A± = B{0,p) f~}(i), A2 = B(x,p)f](ii), 
and A3 = R \{Ai\^jA2). It is straightforward to check that when y G A\, \y — x\ > p/2, so 
J Al ^^7^ £ Io rd ~ l ~ a d r ^ ^ ^J=3- Similarly, we have J^ ^^^dy < ^=j. For 
y G A 3 , we check that \y - x\ > \y\/2, so J M j^j^jdy < f™ ra+ ^ d+1 dr < pa ^_ d . 



For the second inequality, we only have to note that 



f 1 1 f 1 f 1 1 

/ rr^ 1A i n?^ / n^ d y + n^i ^< c 

jRrf lz/r k-yr 7id<2 |yr 7id>2 |yr k-yr 



for |x| < 1. the proof is complete. D 



(5.8) 



Lemma 5.4. For (3 > 0, 



1 , 1 , \y\ 2 



( ll/r (2vrt)2 2t 



as t — > 00. 



1 \ l^ Q ™f ll/l 



2 



Proof. We write / Rd (l A t-™) 3- exp(--!fL)dy = (z) + (z'z), where 

(z) = /" — ^ exp(-^L)dy -> (5.10) 

J\y\<i (2vrt)2 2i 

as £ — > 00, and 



/" 1 1 Iwl /"°° 1 r 

**)=/ U^T ^rexp(-— -)rfy= / r _1_/3 - - r exp(- — )a> 

J\y\>i \y\ P (27rt)3 2t -/l (2vrt)2 2t 

1 d-2-,3 1 Z" 00 d-2-ff 

= j 2 2 — / A 2 exp(— A)aA. 

(2vr)2 t § 4 



(5.11) 



So if ^=§=£ > -1, i.e., < d, {ii) < fl If /3 = d, {ii) < r^ Int. If > d, {ii) < H from 
integration by parts. The proof is complete. □ 

Lemma 5.5. When a G (0, 1), f R2 q £ {x)q £ {y) ^ z+ ^_ y ^ dxdy -)■ ^ as e ->■ /or z / 0. 

Proof. By change of variables, we write 

1 /" 1 

qe{x)q £ (y)- — ■ -dxdy= / q £ {w + y - z)q £ {y)-—dydw 

\z + x-y\ a J R 2 \w\ a 



/ ,+ / ) Qeiw + y - z)q £ {y)——dydw 

J\ W \<\A J\ W \>\A \w\ a 



1 A A ( 5 - 12 ) 

'M<f ^| W |>Hl 



and since 

(ii) = 

'|v^+z|>Jfi 
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(H) = l Q{w + y)q{y)—= — ■ — —dydw, (5.13) 

JWew+z\>& \y/ew + z\ a 



by the Lebesgue dominated convergence theorem, we have (ii) — > j^ as e — > 0. For (i), we write 

00= / ,+ / l )qe{w + y-z)q £ {y)-—dydw. (5.14) 

\J\ W \<^,\ V \>^ J\M<^,\y\< l 4J H a 

For the first term, use g £ (|z|/4) to bound q £ (y), then integrate in y,w; for the second term, use 
g £ (|z|/4) to bound q £ (w + y — z), then integrate in y,w. Since q £ (\z\/4) — > as e — > 0, we have 
(i) — )• 0. The proof is complete. □ 

Lemma 5.6. Assume a G (0, 1), then f R2 QeA^i + yi)q £2 (x2 + y 2 )\yi -y2\~ a dy\dy 2 < C\x\ - x 2 \~ a 
for some uniform constant C. 

Proof. See Lemma A. 2. in [10]. □ 
Lemma 5.7. When d > 3 and a G (0, 2), 

lim — / / P(|J3 a - B tt | < e)dsdu = 0. (5.15) 

Proof. By explicit calculation, we have 

— / / H\B S - B u \ < e)dsdu 
£ Jo Jo 

j ft r foo d ^ 

Az~ 2 e~ A , , , n d\dxds 



(7r)ie« io 4|< £ ' 



r« 7fl ./irl^./Ml X 

2s 

fOO /• /"OO 



/■OO /* /'OO 1 

III 1 kl< £ l|xp<2Asls<tA 5 ~ 2 e~ A _ 2 dActeds ( 5 16 ) 

JO JR d Jo Fl 



(7r)2 £ « Jo 



dA / A^~ 2 e~ A ( AsL e2 + -e 2 F e 2 I l s<t ds 
(ir)is<*Jo J V A< ^ 2 



dAA2~^e _A 1.2 H 1,2 -1,2 = (i) + (u) + (Hi). 

(7r)l £ « 7o V 2 tx>* 2 fx<* 8A fx<*/ 

We check that (i) ~ e d_a , and (ii) ~ e 2_a , (m) ~ e A ~ a + e rf_Q , so the proof is complete. □ 
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